Kondo screening and Magnetism at Interfaces 
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The nature of magnetic order and transport properties near surfaces is a topic of great current 
interest. Here we model metal- insulator interfaces with a multi- layer system governed by a tight- 
binding Hamiltonian in which the interaction is non-zero on one set of adjacent planes and zero 
on another. As the interface hybridization is tuned, magnetic and metallic properties undergo an 
evolution that reflects the competition between anti-ferromagnetism and (Kondo) singlet formation 
in a scenario similar to that occurring in heavy-fermion materials. For a few-layer system at in- 
termediate hybridization, a Kondo insulating phase results where magnetic order and conductivity 
are suppressed in all layers. As more insulating layers are added, magnetic order is restored in all 
correlated layers except that at the interface. Residual signs of Kondo physics are however evident 
in the bulk as a substantial reduction of the order parameter in the 2-3 layers immediately adjacent 
to the interfacial one. We find no signature of long range magnetic order in the metallic layers. 

PACS numbers: 71.10.Fd, 71.30.+h, 02.70.Uu 
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Sufficiently strong electronic correlations can cause the 
formation of an insulating phase at commensurate fill- 
ings. In general, there is a non-zero critical interaction 
strength required for this "Mott transition" , so that two 
uncoupled bands with different degrees of correlation can 
coexist in metallic and insulating states. The behavior 
of spectral functions and magnetic and superconducting 
correlations, when interband hopping or interactions are 
turned on, is a challenging theoretical problem. The cou- 
pling could immediately force both bands to be in the 
same (metallic or insulating) phase, or coexistence might 
persist up to some critical degree of coupling [H 0] . 

Closely related questions arise as clean interfaces be- 
tween correlated materials become accessible Q. Here 
the role of different orbitals is played by the multiple 
layers. It has been suggested that it might be possible to 
"engineer" specific forms of spectral functions at the in- 
terface by varying the materials partnered, as well as de- 
sign other properties arising from electronic interactions 
[flHff . Experimental realizations include tunable 2D elec- 
tron gases in oxide (SrTiOa/LaAlOs) heterostructures, 
control of magnetoresistance at manganite interfaces Q , 
novel magnetic properties at boundaries between cuprate 
superconductors 0] , and observation of magnetic proxim- 
ity effect in Cu/CuO interfaces @. 

While the detailed chemistry of both multi-orbital and 
layered materials is complex, an interesting starting point 
for studying the qualitative properties of metal-insulator 
interfaces is provided by the multi-layer Hubbard Hamil- 
tonian. In this model, electrons have both intralayer and 
interlayer hopping, as well as layer-dependent contact in- 
teractions. The parameter space is large and in this paper 
we focus on the simplest realization of the physics of a 
metal-insulator interface in which all hybridizations are 
chosen to be equal except the one at the interface; the 



corresponding Hamiltonian is 
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Here c\ la {c ila ) are creation (destruction) operators for 
fermions of spin a at site i in layer Each layer is 
an A-site square lattice with a contact interaction Ui 
chosen to be non zero, Ui = U, on "correlated" layers 
I = 1, 2, 3, ■ • -, and zero for an additional set of "metal- 
lic" layers (I = — 1, — 2,— 3, ■ ■ ■). Layers are arranged in 
order of increasing I so that I = ±1 label the layers at the 
interface, t and tw are the intra and interlayer nearest- 
neighbor hybridizations, tw = t except at the interface 
where it takes the value = V. We consider the case 
where \i = which, as a consequence of particle- hole sym- 
metry, makes all layers half- filled, (nu a ) =0.5. Recent 
studies on similar models have found induced magnetic 
order in the metal (lol | and quasi-particle penetration in 
paramagnetic Mott insulators [11 



Questions that arise in connection with Hamiltonian 
([1]) can be seen as extensions to those t ypi cally asked in 
the context of heavy- Fermion materials [12J and concern, 
at least at half-filling, the competition of magnetic order 
and screening of local moments by conduction electrons. 
Heavy-Fermion materials are modeled by Hamiltonian 
([T]), a bilayer with / = ±1, or by its strong coupling limit, 
the Kondo-Hcisenberg lattice, where charge fluctuations 
on the correlated layers are neglected. The fundamental 
issue we address here is how this competition is affected 
as the two-layer case crosses over to the 3-dimensional 
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bulk-to-bulk interface. 

At small and large interface hybridization, the system 
is adiabatically connected to, respectively, the V = 
and V = oo limits. At small V the system is made up of 
magnetically ordered layers weakly coupled to a metal. 
At large V the central bilayer decouples and leaves the 
external layers either metallic (7 < —2) or insulating, 
with anti- ferromagnetic long range order (/ > +2). 

Our results indicate that, for a four-layer system, two 
interacting and two metallic sheets, and at intermedi- 
ate interfacial hybridization, there exists an intervening 
phase where loss of anti-ferromagnetic order is seen in 
both correlated layers, +1 and +2, despite the fact that 
the latter is not in direct contact with the metal. We 
found that the electronic structure of the metal is also 
profoundly affected and that the overall phase of the 
quad-layer can be characterized as a Kondo insulator. 
This is in contrast to our other finding when the inter- 
action region becomes thicker than the metallic one: no 
loss of magnetic order is found in layers beyond the one 
immediately adjacent to the metal, regardless of the hy- 
bridization strength V. In this case, upon increasing V, 
a direct transition between the small and large V regimes 
results. 

We addressed the physics of Hamiltonian (TP) using de- 
terminant Quantum Monte Carlo (DQMC) 13|, an ex- 
act, finite-T method for solving tight binding Hamilto- 
nians on finite lattices. As we limit our calculations to 
the perfectly half-filled case, there is no sign problem at 
any temperature. Our results are averaged over several 
independent simulations, and the error bars correspond 
to the standard deviation of the mean. The imaginary- 
time step is set to At = t/8. We present results for the 
in-plane anti-ferromagnetic structure factor, 
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and the local layer dependent spectral function Ai(ui), 
obtained by inverting the integral equation 
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via the maximum entropy method [14| . Gi(t) = 
Ylitr(T c U<t(t) c !z<x(0)) is the quantity directly obtainable 
by DQMC, and is averaged over 4 boundary conditions 
corresponding to setting the hopping at the boundary of 
each layer to ±i 15|. We also study the in-layer electri- 
cal conductivity ai which is extracted from the current- 
current correlation function 
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with j x (r,0) = itJ2A c t+xM C , 
on the intralayer contribution to A 2 



Wc focused 
only, assuming 



that this correctly characterizes the conductive property 
of each layer, oi is extracted using the approximate form 
of the fluctuation-dissipation relation, valid at large /3 
and first discussed in I16II. 



A ra , z (fc = 0,r = /3/2) = nai/(3 2 



(5) 




FIG. 1: (color online) (a) V-dependence of local moments, 
mi = ^2 (( a ii) 2 ) /N, on each layer, when two metallic layers 
are coupled to two correlated ones, (b) Finite size scaling 
of in-plane structure factor S| of the correlated layer (I = 
2) farthest from the interface. For small V, there is long 
range order in the thermodynamic limit which vanishes for 
intermediate V and is recovered for large V. S% reaches its 
ground state value at fit = 10. 

To gain initial quantitative understanding of the evo- 
lution of magnetic properties as V increases, we show, 
in the top panel of Fig. [TJ the evolution of local mo- 
ments in a system of two metallic and two interacting 
sheets. There are three regimes, most clearly evidenced 
by the behavior of the metallic layer at the interface. At 
V < t local moments on layer —1 are essentially identi- 
cal to those of a non-interacting system. In t < V < 4t 
the moments monotonically increase and they saturate 
at V ~ At. The evolution of the other layers follows nat- 
urally, with layer +1 merging with —1 at large V in a 
phase that can be best characterized as a band insula- 
tor made of weakly interacting dimcrs. Layer +2 has the 
only non-monotonic evolution: magnetism is first sup- 
pressed and then revived as the central dimer phase gets 
increasingly stabilized. 

We then use finite size scaling on 5f f to investigate 
whether order in layer +2 is lost in the regime where 
the moments are most suppressed. As shown in the 
lower panel of Fig. [TJ 5| scales to a nonzero value for 
both small and large interface hopping, V, when plotted 
against the inverse linear system size 1/y/N. For both 
these regimes, there is long range order in the ground 
state in the thermodynamic limit 17] as expected from 
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FIG. 2: (color online) (a) Spectral function A2(u)) in corre- 
lated layer I — 2. The Slater gap present at small V due to 
AF order vanishes at intermediate V and reappears at larger 
V. (b) In correlated layer 1 = 1, Ai(ui) resembles ^2(01) but 
unlike I = 2 the large V behavior is a broader gap associated 
with the singlet energy scale. In both non-interacting layers 
A-i(u>) (c) and A-2{uf) (d) a gap opens as V increases. At 
large V, layer —2 recovers metallic properties while the singlet 
gap is visible in layer — 1. 




FIG. 3: (color online) (a) In-plane conductivity <r_i in the 
metallic layer / = — 1 as a function of V, at several inverse 
temperature ft values. At intermediate V is small but non- 
zero, but vanishes at V > At due to dimer formation, (b) In- 
plane conductivity cr„2. This noninteracting layer becomes 
insulating for intermediate hybridization, t < V < 3t and 
then recovers when the pairs are fully formed and pinned at 
the interface. 



the behavior of the local moments. However, in the inter- 
mediate regime (starting at V > t, in rough correspon- 
dence to where loss of anti-ferromagnetic order happens 
in the periodic Anderson model [3] and the bilayer Hub- 
bard model ) loss of magnetic order is clearly observed 
on layer +2 as well. 

The most likely candidate mechanism for such loss of 
order involves Kondo screening in both interacting lay- 
ers. As the formation of a resonance in the single particle 
spectral density is one of the hallmark of such process, 
we plot, in FigFJJ the layer dependent spectral density 
at T = t/30. We found that, at V = 2t, both interact- 
ing layers are characterized by the presence of a Kondo 
resonance. Due to the fact that we are focusing on a half- 
filled system, the resonance is split as typically happens 
for Kondo insulators. 

We can gain further insight into the nature of this in- 
termediate phase by looking at the behavior of the non- 
interacting layers. Figure [3] shows the conductivity in 
layers I = — 1 and I = — 2. At small hybridization, before 
the loss of magnetic order, the conductivity increases as 
T is lowered, showing these two sheets to be metallic. 
In the intermediate regime t < V < 3i, the conductiv- 
ity is strongly suppressed in both layers, and our inverse 
temperature results suggest that these layers become in- 
sulating around V = t in correspondence with the loss of 
magnetic order. 

From these results we can draw a few significant con- 
clusions. First, that there is a Kondo proximity effect, 



as already observed using dynamical mean-field theory 
on a similar model (llj . The novelty of our finding re- 
sides in our treatment of non-local correlation. It allows 
for a proper description of the competition between mag- 
netic order and Kondo screening and shows the latter to 
be effective in destroying order even on layers not di- 
rectly coupled to the metal. Furthermore, states from 
both metallic layers participate in the screening of local 
moments as best evidenced by the drop in conductivity. 
This situation is reminiscent of a long-range resonating 
valence-bond state although it is unclear whether such a 
description remains meaningful in the present context of 
itinerant electrons. 

The other kind of proximity effect that is expected in 
such systems is due to the presence of the magnetically 
ordered layers at smaller values of V. This would seem 
a likely scenario, especially since the metallic phase lives 
on layers with nested Fermi surfaces with infinite, T = 0, 
anti-ferromagnetic susceptibility. Our calculations, how- 
ever, do not find any significant penetration of magnetic 
order in the metallic layers. Although we cannot exclude 
that an extremely small order parameter might develop 
at low T for some range of V, our finding suggests that 
such an order would not survive the generic scenario of a 
system with finite susceptibility. This result appears at 
odds with recent experiments [9( finding evidence for an 
anti-ferromagnetic proximity effect. 

We now consider the question of how the competition 
of Kondo screening and magnetic order is affected when 
additional interacting layers I = 3, 4, ■ • ■ arc present. The 
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FIG. 4: (color online) (a) In-plane structure factors Sf { as 
functions of V for a system of four correlated layers coupled 
to two metallic ones. At intermediate V antiferromagnetic 
order is suppressed in every correlated layer. At large V the 
order recovers in every correlated layer but the one at the 
interface, (b) Finite-size scaling of the in-plane structure fac- 
tors Sf for a system of six correlated layers coupled to two 
metallic ones, at V = 2t. A systematic reduction of spin corre- 
lations is evident as the metallic interface is approached from 
the correlated side. Long range magnetic order is completely 
destroyed in the correlated layer at the junction. 



of the metal and strongly interacting material is a sup- 
pression of magnetic order on the correlated side. We 
did not observe the converse phenomenon, namely a sig- 
nificant penetration of magnetism into the metal, as has 
been noted in ■ It is possible this difference arises from 
the lower value of the on-site interaction, U/t = 4, used 
here, compared to U/t w 17 in [Io[. Such large couplings 
are difficult to treat in DQMC. The on-site interactions 
studied here are relevant to the range of fitted values for 
a number of interesting materials, e.g. CuO 
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We showed that for thin insulating layers the Kondo 
effect embraces correlated and metallic layers that are 
not in direct contact with each other. Although such an 
extended Kondo insulating phase does not require any 
fine tuning, it is not a phase that permeates a large frac- 
tion of the "phase space" for such systems. Here we have 
shown, for instance, that forming Kondo singlets across 
multiple layers is a process that can be defeated by mag- 
netic proximity effect of layers farther from the interface. 

This work was supported by: the CNRS-UC Davis 
EPOCAL LIA joint research grant; by NSF-PIF-1005503 
and DOE SSAAP DE-FG52-09NA29464; ARO Award 
W911NF0710576 with funds from the DARPA OLE Pro- 
gram; and Army Research Office Grant 56693-PH. 



resulting system can be thought of describing the in- 
terface between a thin metallic film and a bulk anti- 
ferromagnetic Mott insulator or as a heterostructure 
where the insulating domains are substantially thicker 
than the metallic ones. A scan of 5f f for different values 
of V on a 6-laycr cluster with 600 sites (top panel in Fig. 
2]) indicates that the most likely parameter regime to ob- 
serve loss of magnetic order in layer 2 is for V/t G [1,2]. 

However, for V = 2t (lower panel of Fig. 2]), we found 
that, while Sf does not extrapolate to a non-zero value, 
the same does not happen for layers located deeper into 
the interacting material i.e. I = 2,3,4,5,6 are all anti- 
ferromagnetically ordered. This revival of magnetic long- 
range order on layer +2 can be interpreted as a magnetic 
proximity effect exerted by a bulk anti-ferromagnet on 
those correlated layers subject to Kondo screening. Note 
that a clear suppression of the order parameter is still 
observed for I = 2 and several layers deeper in indicat- 
ing coexistence of Kondo screening and magnetism (the 
anomalously large value at I = 6 is a known phenomenon 
where surface magnetic correlations are larger than the 
bulk and has been widely explored experimentally [20l]). 

In conclusion, we have presented results on a model of 
metal-insulator interface, the multilayer Hubbard Hamil- 
tonian. Even within this simplified tight-binding model, 
there are many possible choices of the intralayer and in- 
tcrlayer hoppings. The dominant feature of the coupling 
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